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Abstract 

The chiral Gross-Neveu model or equivalently the linearized Bogoliubov-de Gennes equation has been mapped to the nonlinear 
Schrodinger (NLS) hierarchy in the Ablowitz-Kaup-Newell-Segur formalism by Correa, Dunne and Plyushchay. We derive the 
general expression for exact fermionic solutions for all gap functions in the arbitrary order of the NLS hierarchy. We also find 
that the energy spectrum of the «-th NLS hierarchy generally has n + 1 gaps. As an illustration, we present the self-consistent 
two-complex-kink solution with four real parameters and two fermion bound states. The two kinks can be placed at any position 
and have phase shifts. When the two kinks are well separated, the fermion bound states are localized around each kink in most 
parameter region. When two kinks with phase shifts close to each other are placed at distance as short as possible, the both fermion 
bound states have two peaks at the two kinks, i.e., the derealization of the bound states occurs. 

Keywords: Bogoliubov-de Gennes equation, Chiral Gross-Neveu model, Nonlinear Schrodinger hierarchy, AKNS formalism 
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1. Introduction 

The Bogoliubov-de Gennes approach jlj] was originally for- 
mulated for dealing with spatially inhomogeneous supercon- 
ductors, that consists of solving the Bogoliubov-de Gennes 
(BdG) equation and the gap equation self-consistently. It has 
been widely applied from condensed matter to high-energy 
physics to describe various phenomena such as solitons in 
charge-density- wave systems J2l, conducting polymers 

incommensurate spin-density-wave states in doped Cr 
_ and fermion number fractionization in relativistic quan- 
tum field theory 0, [Toll . In particular, an important application 
in high-energy physics is to the (chiral) Gross-Neveu model 
fllll and equivalently the Nambu-Jona-Lasinio model fl2ll that 
are known to describe dynamical chiral symmetry breaking and 
resulting dynamical mass generation. The original BdG sys- 
tem reduces to these relativistic field theories by linearization 
known as the Andreev approximation in the theory of super- 
conductivity. 

However, it is generally a difficult task to find analytic self- 
consistent solutions of the BdG and gap equations when the 
gap function is inhomogeneous. For real-valued gap functions 
in one-dimensional systems, analytic solutions were completely 
known under uniform boundary conditions at spatial infinities; 
a real kink fl3l l2ll. a kink-anti-kink bound state (polaron) lfl4l 
15], a three-real-kink Jl6l 17 1 and more general real solutions 



flill . Under a periodic boundary condition, a real kink crystal 
known as the Larkin-Ovchinnikov state is known lfl9lEoll . On 
the other hand, few analytic solutions in complex condensates 
are known; a twisted (or complex) kink 12111 . i.e., a kink with a 
twisted boundary condition on the phase. 

Recently, a breakthrough has been made by Bagar and Dunne 
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1 22112311 that a suitable ansatz for the Gorkov resolvent reduces 
the gap equation to the nonlinear Schrodinger (NLS) equa- 
tion, which is exactly soluble. Since the derived NLS equa- 
tion is a closed equation for the order parameter A(x), this ap- 
proach enables one to avoid the self-consistent calculation of 
the coupled BdG and gap equations. They also found a new 
exact self-consistent crystalline condensate, a complex kink 
crystal with both phase and amplitude modulating (the Larkin- 
Ovchinnikov-Fulde-Ferrel state Jl9ll24ll ). that includes all pre- 
viously known solutions as special cases. This approach by 
Basar and Dunne has been extended to incorporate spin imbal- 
ance effect which plays an important role in superconductors 
under a strong magnetic field and spin polarized ultracold Fermi 
gases 125(1 . 

The Basar-Dunne approach has been further developed by 
Correa, Dunne, and Plyushchay [2fi] by introducing the inte- 
grable nonlinear equations for the gap function A(x) that be- 
long to the NLS hierarchy of the celebrated Ablowitz-Kaup- 
Newell-Segur (AKNS) formalism. However, general expres- 
sion for fermionic solutions of the BdG equation has not been 
obtained for those A(x); One had to solve the BdG equation to 
obtain fermionic solutions for each given order parameter A(x). 

The aim of the present Letter is to derive the general expres- 
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sion for exact fermionic solutions of the BdG equation when 
the gap function A(x) obeys the NLS equations in the arbitrary 
order of the NLS hierarchy. We first formulate the BdG system 
in the AKNS form and briefly review the derivation of higher 
order NLS equations. We then discuss the derivation of the gen- 
eral expression for fermionic solutions in detail. As a concrete 
example, we apply our formalism to the two-complex-kink so- 
lution for the second nonlinear equation (AKNS2) to derive the 
fermionic solutions for the bound states as well as for scattering 
states. 

2. AKNS form and solutions of linearized BdG equation 

2.1. Problem considered in this Letter 

First, let us clarify the problem treated in this Letter mathe- 
matically. What we want to solve is the linearized (Andreev- 
approximated) BdG equation 
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(2.1) 



with the assumption that the gap A(x) satisfies stationary higher 
order NLS equations B27I1 . known as AKNS,, (n = 1,2,... )|2g|: 



n+2 

iA f = £c/-iM<f), 



(2.2) 



where cjs are arbitrary real coefficients, and M^s are given by 



-iMg = A, 
-iMfl = — iA t , 



-iM 



(31 



-A xx + 2|A| 2 A, 



-iM® = i(A xxx - 6\A\ Z A X ), 



-iM 



(5) 



2(|A| 2 ).„A - 3A*(A 2 )„ + 6|A| 4 A, 



(2.3) 
(2.4) 
(2.5) 
(2.6) 
(2.7) 



with 2x2 matrices U and V defined below, and the nonlinear 
equation that we want to solve arises as a compatibility condi- 
tion (or a zero-curvature condition) 

U t -V x + [U, V] = 0. (2.10) 

In the case of NLS hierarchy, U is a 2 x 2 matrix defined by 

-iA q(x, t)\ 



U{x, t, A) 



r(x, t) 



iA 



(2.11) 



Here A is an (x, ^-independent spectral parameter. Henceforth 
we often omit arguments of functions when it is clear from the 
context. If we identify r, q, and A as 



-iA, 



iA*, and A 



(2.12) 



then Eq. (12.8b readily reduces to the BdG equation ( 12.11 ). The 
form of V for higher order NLS equations will be given later. 

Note that if we consider a time-independent problem, the 
zero-curvature condition has the same form with the Eilen- 
berger equation: 



V X = [U,V]. 



(2.13) 



Therefore, V - R0-3 immediately gives a particular solution for 
the Eilenberger equation up to normalization, where R is the 
Gorkov resolvent. (Note that R must be normalized as detR = 
-l/4§i.) 

Even though V is an ^-dependent matrix, the eigenvalues of 
V is independent of x, since 



(trV), = (trV 2 ), = 



(2.14) 



follows from Eq. (12.13b . and therefore a characteristic polyno- 
mial for y 



det (y - iu>l 2 ) = det V - iu> tr V - of 



Here the subscripts t and x denote the differentiation with re- 
spect to t and x, respectively. Since we consider a stationary 
(time-independent) problem, henceforth L.H.S. of Eq. (12.21 1 is 
set to zero. How to generate A/gs is reviewed in Section [231 
For later convenience, we also prepare the term "pure AKNS," 
for the AKNS,, with c\ — ■ ■ ■ — c n+ \ = and c n+ o - 1. 

The most familiar NLS equation is equivalent to AKNSi. 
AKNS2 is also known as the Hirota equation 12811 and used to 
describe the soliton propagation phenomena in optical fibers. 
(See also Ref. lf29ll .) We also note that pure AKNS2 is equiva- 
lent to the modified KdV equation when A is real-valued. 

2.2. Overdetermined system and Compatibility condition 



The AKNS system H30tl is formulated through the following 



overdetermined system: 

8_lu 
dx \v 



U(x, t, A) 



(2.8) 
(2.9) 



1 2 2 2 ( 2 - 15 > 

= -(tr y 2 - (tr V) 2 ) - ia> tr V - ar 



becomes independent of x. It is worthy to note that Eq. ( 12.131 ) 
is identical to the famous Lax equation |[3~ll] if we replace x by 
t, and the above proof is the same as the proof of the isospectral 
property of the Lax operator. 

2.3. Higher order NLS 

In this subsection we review how to generate higher order 
NLS equations. Although the content in this subsection is only 
a revisit of the famous textbook by Faddeev and Takhtajan 12711 . 
we briefly review it for self-containedness since the definition 
of spectral parameter in our Letter differs from theirs in twice 
factor. 

The NLS equation has infinitely many conservation laws, and 
correspondingly, there are infinitely many integrable equations. 
In accordance with Ref. ll26Tl . we call the equation generated 
from the n-th order conserved quantity AKNS„_2- Let be a 
matrix which yields pure AKNS„_2- Then V (n) can be calculated 
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as follows 1 2711. Let M be a 2 x 2 matrix which satisfies the 
equation 



M x = [U,M] 
and has the following formal Laurent series 

^ M (,,) rm io-3 

M = > , M (0) : = 

^ (-2/1)" 2 



(2.16) 



(2.17) 



(2.18) 



Substituting Eq. I l2.17l l into Eq. ( 12.16b . and using 

U = 2AM^ + U (l \ E/ (1) :=f° 
one obtains the recurrence relation 



M< n) = ^[cr 3 ,M ( " +1) ] + [t/ (1) ,M (,,) ] (n = 0,l,...). (2.19) 

Thus one can determine M (,,) inductively. If one sets all arising 
integration constants to zero, one obtains 



- { -r "oV 



-ir x lrq 



M (3) _ [ Vx ~ rq x q xx ~ 2rq 
\r xx - 2r 2 q -{qr x - rq x ) 



( Krq X x + qr xx 



M' 4 » = 



-r x q x - 3rV) 



-Kq.xxx - 6rqq x ) 



■i (. \ ~Krq xx + qr xx 
<r xxx -6qrr x ) -r x q x -3rW) 



(2.20) 
(2.21) 
(2.22) 

, (2.23) 



y(") f or 

pure AKNSn-2 is then given by 



= J^(-2A)"- l - k M (k) . 



(2.24) 



k=0 



Setting V = V (n) in zero-curvature condition ( 12.101 ) and using 
( 12.19I ). one obtains the higher order NLS equation 



U) 1 ' = ~[cr 3 ,M m ) ^ 



(») 



An) 



(2.25) 



Thus, the off-diagonal element of M (/,) gives the higher order 
NLS equation itself. The matrix V which yields general (non- 
pure) AKNS,, is simply the real-coefficient linear combination 

of V (J) s: 



n+2 



(2.26) 



and the resultant equation is given by iq t = 2y=i c j(~^n^ as 
shown in Eq. (12.21 1 with recalling q = -iA. 



2.4. Solutions for stationary AKNS equation 

In this subsection, we solve Eq. ( 12.81 ) in static case, namely 
the ordinary differential equation 



(2.27) 



under the assumption that there exists another matrix V(x, A) 
satisfying the differential equation 



V X = [U,V]. 



(2.28) 



Although our main interest in this Letter is the case where U 
and V are given by Eqs. d2.1 II ) and ( 12.26b . the solution given in 
this subsection does not depend on the special form of U and 
V. 

The solution is constructed from the following ansatz: 



(2.29) 



The ansatz d2.29b and the original equations ( 12.271 ) and ( 12.28b 
are "compatible", since ia> x (Eq. ( 12.27b ) and <9, t (Eq. ( 12.29b ) 
yield 



respectively, while R.H.S. of these two are equal because of 
Eq. (12381 . 

Even though we can treat the problem in the above frame- 
work for general U and V, we restrict our discussion to the case 
where U and V are traceless: 



tr U = tr V = 0, 



(2.31) 



since the final expression for (if, v) becomes a bit simpler. We 
note that the matrices U and V for the NLS hierarchy, Eqs. 
( 12. 1 lb and ( 12.26b satisfy this condition. 

For later reference, we write down the time-independent 
zero-curvature condition ( 12.28b for each component: 

Vux = -V 22x = U n Vn - U 2i V n , (2.32) 
V nx = 2t/iiVi 2 - 2U n V n , (2.33) 
V 2 i x = -2UuV 2 i+2U 2l Vu. (2.34) 

Here the traceless assumption is used. 

In order for Eq. (12.29b to have a non-vanishing solution, 

det(V - i(ol 2 ) = <-> a? = det V = -Vr, - V l2 V 2 \ (2.35) 



must hold. As discussed in Section l2T2l det V is independent of 
x, and u> is determined as an x-independent constant: 



CO 



= ± Vdet V. 



(2.36) 



From the expression (12.35b we also obtain an important factor- 
ization: 



Vi 2 V 2 i = (iV u + a>)QVu - a). 



(2.37) 
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This relation technically plays an important role in the next sub- 
section. When co is a root of Eq. J2.35I ). 



iV 



21 



(2.38) 

// IV12 IV 11 - co 

holds from Eq. ( 12.291 ). Using this, one can eliminate either of 
u or v. If one chooses to eliminate v, from the first row of Eq. 
( 12.27b . one obtains a first order differential equation for u : 



u x — U\\u + UylV = 



iVu + co 



iV 



12 



With the use of Eqs. ( 12.331 ) and ( 12.37b . one can show 

2icoUn V*2i 

u 



2 «£ _ Vnx 

V n (iVu + co)(iV u - co) ' 



Furthermore, since it follows that 

TJ T/ U 12 V 2 i + U 21 Vi2 Vn,- 
UnVn = + -7— 



(2.39) 



(2.40) 



(2.41) 



from Eq. ( 12.32b . one obtains 

itoVux 



2 ± = Vnx + 



Uu U 2i 

+ ICO I + 

Vn (iVu + co)(iV n - to) \Vn V 2 \ 



(2.42) 



Except for the rightmost term, each term can be integrated sym- 
bolically, and the solution is given by 



, \\V\l - CO 

vr = CV n A / — exp 

1 Vl I + CO 



r , (Uu u 2l 

i dx 1 

Jo \ Vn V 21 



(2.43) 



and using Eq. ( 12.38b . the expression for v is also obtained: 



iVu+aj 

= -CV 2i A — exp 

lVl i — co 



ico 



Jo l Vn V a 



(2.44) 



Here C is an integration constant. Note that there are two possi- 
ble values for co since Eq. ( 12.35b is quadratic, and it corresponds 
to two linearly-independent solutions of the original equation 
(12.27b . (The case co — is rather exceptional and the second 
solution must be constructed by reduction of order.) Needless 
to say, when one takes the square root of (12.43b and ( 12.44b . one 
must choose the sign so that the relation ( 12.38b is satisfied. 

The solutions ( 12.43b and ( 12.44b are one of the main results 
of this Letter. Since this expression is valid for arbitrary spec- 
tral parameter A, it means that the linearized BdG equation has 
been solved for arbitrary energy E regardless of whether the 
wave function of the energy E diverges or not. In the next sub- 
section, we give a criterion for a given energy E to belong to 
the spectrum. 

2.5. Energy spectra 

In the preceding subsection, we have derived the general so- 
lution for the stationary AKNS equation. Here we concentrate 
on the specialized issues for the linearized BdG equation, i.e., 
we consider the case where the matrices U and V are given by 
Eqs. (f27TTT ) and d2?26t . 

Henceforth we always assume q = -iA and r = iA*, so 
r = q*. As seen in Section [231 V of NLS hierarchy always 
satisfies: 



(i) If A is real ( <-> E — -A is real), iVu is real. 

(ii) If A is real, V*, = V n . 



From these facts, the left hand side of d2.37b is real and non- 
negative when A is real. We can further see the following: 



(iii) If A is real, det V is also real. Therefore, co — ± Vdet V is 
either real or pure imaginary. 

(iv) If A and co are real, \V\\ + co and W\\ - co are both real 
and have the same sign. Since iVu - \co\ < iVu < iV\\ + 
\co\, W\\ also has the same sign. Furthermore, if co + 
and the gap function A(x) is bounded and C k with k being 
sufficiently large, sgn(iVn) is globally constant and does 
not depend on x. 

(v) As a corollary of the above, if A and co are real, the ratio 
(iVu - co)l(W\\ + co) is real and non-negative. 

We need an explanation for the latter part of (iv). If iV\\ + \co\ 
changes its sign at some point x, then the sign of iVu - \co\ also 
must change simultaneously. However, it is impossible since 
iVu is continuous and \co\ > 0. 

From the above information (i)-(v), let us estimate the 
asymptotic behavior of Eqs. ( 12.43b and (12.44b . When A = -E 
is real, the integrand in the exp term ky 1 + tt- = 2 Re ^ is a 
real-valued function. Therefore, if co is real, the absolute value 
of the exponential term is equal to unity. Furthermore, using 
the fact (v) and Eq. (12.37b . when A and co are real, the relations 



M = |C||iVi 
follow, because, e.g., 



IV 



121 



iVu - co Wu-co 

= A/V12V21— = |iVn -co\ 



(2.45) 



(2.46) 



iVu + co \ iVn + co 

holds. The expression (12.45b is obviously bounded if A(x) is 
a bounded function. Thus, if co — + Vdet V is real, the corre- 
sponding energy is in the spectrum. On the other hand, if A is 
real but co is pure imaginary, the exp term diverges exponen- 
tially. To sum up, 

(A) If co 2 = detV > for a given E, this energy is in the 
spectrum. 

(B) The band edges are determined by the equation co 2 = 
det V = 0, which is an equation of order 2n + 2 with re- 
spect to A (or E) if the gap A(x) is a solution of AKNS„. 
So there exist 2n + 2 band edges in general. In this case, 
the solutions ( 12.43b and ( 12.44b are drastically simplified: 



("V) 



(Vi2,-V 2 i), 



(co = 0). 



(2.47) 



(C) If co 2 = det V < for a given E, this energy is not in 
the spectrum. The solution diverges exponentially and un- 
physical. 

Since the discrete spectrum can be regarded as a shrinking limit 
of continuous spectrum, as a special case of (B) we arrive at: 
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(D) If the equation det V — has a multiple root, this root be- 
longs to the discrete spectrum. The corresponding eigen- 
state can be a normalizable localized mode. 

Here we give a few technical remarks. If the gap function 
A(x) is a solution of AKNS,,, it also can be a solution of higher 
order AKNS m (m > n). However, when we calculate the posi- 
tion of band edges, we should use the matrix V for the lowest 
order AKNS,,. If we use a higher order one, we might obtain 
additional dummy edges and dummy discrete eigenvalues. 

If A and a> are real, in other words, if E = -A is in the spec- 
trum, one can explicitly write down the square root of ( 12.43b 
and (12.44b with the use of Eq. ( 12.45b . It is given by 



u 

vl ~ \isgn(iV n ) VjiVn + oj 
ho 

y 



are V, 



x exp 



Vl2 V 2 i I 



(2.48) 



where one should recall (iv), i.e.. 



sgn(iVn) = sgn(iVn + oj) - sgn(iV n - oj) 



(2.49) 



is independent of x and this sign factor is necessary to satisfy 
the relation (12.38b . 

3. Example: Complex Two-kink solution in AKNS2 

In this section, as an example, we consider the complex- 
valued two-kink state in AKNS2 and explicitly write down 
fermionic wave functions (u, v) in this state. Let V be 



V = c I V (1, + C2 V (2, + c 3 V (3) + V (4 \ 



(3-D 



where we set C4 = 1 without loss of generality. We then obtain 
the fundamental equation 

ciq - ic 2 q x + c 3 (-fc + 2\q\ 2 q) + i(q xxx - 6\q\ 2 q x ) = (3.2) 



with setting r — q*. We can construct the following three inde- 
pendent integration constants by the method given in Ref . 113211 : 



Ji = YiJijCj, (*= 1,2,3), 



(3.3) 



7=1 



(3.4) 



where J,j is a skew-symmetric 4x4 matrix with components 

Jn = rq, J i3 = i(r x q - rq x ), 
Ju = 3r 2 q 2 + r x q x - r xx q - rq xx , 
J23 = -r 2 q 2 + r x q x , J 1A = i(r xx q x - r x q xx ), 
Ju = 4rV + r 2 q 2 x + r\q 2 - 1rq{rq xx + r x q x + r xx q) + r xx q xx 

Using J 1, J2, and J3, we can show the constancy of det V: 

det V =16/ - 16c 3 ^ 5 + 4(2c 2 + c\)A 4 - 4(ci + c 2 c 3 )A 3 
+ (c\ + 2c y ci- AJ { )A 2 + (2J 2 + 2c 3 Ji - c { c 2 )A 



■ c 2 7i - C3J2 - Ji- 



As a special case, we consider the asymptotically uniform 
boundary condition 



iq 



m (x — > -00), 

me" M (x -> +00). 



(3.6) 



This boundary condition gives a complex-valued two-kink so- 
lution. In this case, we obtain 



ci = —2c 3 m 2 , J\ — m 2 C2 + 3m 4 , 
J 2 — m 4 C2, J 3 — m 4 C2 + 4m 6 , 
det V = (A 2 - m 2 )(4A 2 - 2c 3 A + c 2 + 2m 2 ) 2 . 



(3.7) 



Introducing 0\ and 62 by the relations 

c 2 = 2;72 2 (2 cos#i cos# 2 -l), c 3 = -2m(cosf?i +cos0 2 ), (3.8) 
with 9 — 6\ +62, then det V is factorized as 

(3.9) 



det V = 16(A 2 - m 2 )(A + m cos 6i) 2 (A + m cos 6» 2 ) 2 . 



Therefore the two-kink solution has a continuous spectrum for 
\E\ > m and two bound states at E = m cos 6\ and m cos 62. 

The two-kink solution can be written as follows. Without 
loss of generality we assume < Oj < n (j = 1,2). Introducing 
the notations 



K\—m sin 6\ , k 2 = m sin 62, a — 
h(x) = 
h 2 (x) = 



2 -\JkiK2 



im(e 161 
-aci(1 + e" 2 ^^"^) + ot '\fic\K 2 

(1 + e -2Kdx~x t ))(l + _ | a |2 ' 



(3.10) 



(1 + e 



-2k\ (x-xi) 



)(1 + e-2^(A-. t2 )) _ \ a \2' 



with real constants x\ and X2 related to kink positions shown 
below, A(x) can be expressed a£] 



A{x) = m + 2i(e^ iSl /ii(x) + e~ i62 h 2 (x)). 



(3.11) 



The normalized bound state for E = m cos 8j (j — 1 , 2) is 
given by 



Uj(x) 
Vj(x) 



jit] \eP'hj{xy) 
and the scattering state with eigenenergy E is given by 



(3.12) 



"(*. E )\ = Jkx 

v(x,E) 



E-kj ' " me 10 ' 
7=1.2 



hj(x) 
E-k \j e 'hj(x)* 



k=± V£ 2 - m 2 . 



(3.13) 



We note that any two-kink solution can become self-consistent 
when the number of flavors is sufficiently large lf33h . 



1 The definition of Xj (i = 1, 2) in this Letter differs from that in Ref. |33 
Here, we use = y [Kje. K ^ x ~ x t ) instead of ej(x) = & c fo~ x fl for |33|| . 



The two-kink solution given by Eq. (13.1 lb has four parame- 
ters, &i,02,x\, and X2- If we exclude a trivial translational de- 
gree of freedom, there are actually three parameters, that is, 
6\ , 2 , and #2 - x\. The two kinks are well separated from each 
other if the following condition holds 

m|x 2 - xi I » 1 or |0i-0 2 |<k1. (3.14) 

In this case, an approximate position of each kink is given by 



Xleft - 



bright — 



log(e 



2 ) 



(17 m 2 -&k 2 )5 2 



32a-- 1 cosh" k(x2~xi ) 

(x 2 »x,) 
(x 2 xi). 



(x 2 » x x ) 

(X2 - Xi ), 



(3.15) 



I z 2K2 to l-costf 

|ii og |^ + '°^; 1 ^) 

I, 2k to l-cos(5 2k 



(3.16) 



Here, we have defined 5 — 0\ - 62 and k = n?sin(0/2). The 
gap function satisfies A(x) ~ mer 2l6> for x\ e n «: x <K jtnght and 
A(je) ~ me~ 2i(fll+fl2) = me -2 '" for x rig ht «jso that -20,- can be 
identified with the phase shift of the 2-th kink. 

Figure Q] shows the complex-two-kink solutions and their 
bound states for various cases. When the two kinks are well 
separated, the fermion bound states are localized around each 
kink in most parameter region. When two kinks have phase 
shifts close to each other, distance between them is shortest for 
x\ - X2- Only around this parameter region, i.e., only when 



m\x2 - X\\ < 1 and \0\ - 2 \ <sc 1 



(3.17) 



hold, the both fermion bound states have two peaks at the two 
kinks, that is, the "derealization" of the bound states occurs. 
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Figure 1: (Color online) Two-kink solution of the gap |A(.v)| 2 (black solid lines) 
and fermionic bound state wave functions, 2(|«i(jc)| 2 + K'i(jc)| 2 ) (red dashed 
lines) and 2(|M2(x)| 2 + |v'2U')l 2 ) (blue dotted lines), for various parameters (we 
set m = 1). The bound states are delocalized around the two kinks when 
phase shifts of kinks are close to each other and xi = x\ [(a)]. As increas- 
ing |*2 - vil with fixed \6\ - 6i\, the distance between the kinks increases and 
consequently the bound states are localized on each kink [(a)— »(b)— »(c)], As in- 
creasing |6i - &2\ with fixed \xj - jci|, the distance between the kinks decreases 
while the derealization character of the bound states is preserved [(a)— »(d)]. 
In (a)— >(e)— >(f), the kink profile is almost unchanged but the derealization of 
bound states is lifted. Thus, the derealization phenomenon of bound states is 
sensitive to the choice of kink-parameters. 



4. Summary and Discussion 

We have analytically derived the general expression for ex- 
act fermionic solutions (12.43t and (12.44b to the BdG equation 
for all gap functions A(x) in the arbitrary order of the NLS hi- 
erarchy in the AKNS formalism, originally suggested by Cor- 
rea, Dunne and Plyushchay lf26ll . Depending on the sign of 
co 2 = det V, these solutions are inside an energy band in the 
energy spectrum (det V > 0), on a band edge (det V = 0), or 
exponentially divergent and unphysical (det V < 0). The en- 
ergy spectrum of our fermion solutions for the AKNS„ contains 
two continuums for scattering states and n bands for localized 
states. Consequently, there exist 2n + 2 band edges, on which 
the fermion solutions become rather simple as in Eq. ( 12.471 ). 
When det V = has a multiple root, the corresponding fermion 
solution belongs to the discrete spectrum. 

As an illustration of our formalism, we have presented the 
complex-two-kink solution for the AKNS2 with the uniform 
boundary conditions for A(x) at x — > ±00. The two kinks 
can be placed at any position and have phase shifts. In this 
example, the equation det V = has two multiple roots, giv- 
ing two fermions localized on each or both of the two kinks. 
When the two kinks are well separated, the fermion bound 
states are localized around each kink in most parameter region. 
When two kinks have phase shifts close to each other and are 
placed at distance as short as possible, the both fermion bound 



states have two peaks at the two kinks. The self-consistency of 
the solutions is shown in 1I33I1 in which the «-kink solution is 
given. While the uniform boundary condition is assumed for 
this example, which can be studied by the inverse scattering 
method, let us stress that our fermionic solutions can be applied 
to any solution under arbitrary boundary conditions. We may 
extend our formalism to unconventional superconductors such 
as multi-gap superconductors and p-wave superconductors. 
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